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Crawford F. Smith and Steve D. Podleski 
Sverdrup Technology, Inc. 

Lewis Research Center Group 
Brook Park, Ohio 44142 


Summary 

The proper use of a computational fluid dynamics code requires a good 
understanding of the particular code being applied. In this report the application 
of CFL3D, a thin-layer Navier-Stokes code, is compared with the results obtained 
from PARC 3D, a full Navier-Stokes code. In order to gain an understanding of 
the use of this code, a simple problem was chosen in which several key features 
of the code could be exercised. The problem chosen is a cone in supersonic 
flow at an angle of attack. The issues of grid resolution, grid blocking, and 
multigridding with CFL3D are explored. The use of multigridding resulted in a 
significant reduction in the computational time required to solve the problem. 
Solutions obtained are compared with the results using the full Navier-Stokes 
equations solver PARC3D. The results obtained with the CFL3D code compared 
well with the PARC3D solutions. 
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1.0 Introduction 


The analysis of the flow in an aircraft inlet, such as the F-18 inlet, at 
subsonic speeds and high angles of attack requires the inclusion of the external 
flow about the forebody, Leading Edge Extension (LEX), and wing in order to 
account for upstream disturbances such as flow separation and shed vortices 
which might be entrained by the inlet flow. 

The numerical solution of this problem is very difficult and requires large 
amounts of computational time. For adequate geometry resolution, grid block- 
ing is necessary. The use of multigridding can sometimes significantly decrease 
the amount of computational time required to obtain a converged solution. In 
addition, proper grid resolution is needed to capture the details of a very complex 
flow field. 

A computational fluid dynamics (CFD) code that has been used to address 
the problem of forebodies at high angles of attack is the CFL3D code [Ref. 
1]. This code has been developed at the NASA Langley Research Center and 
solves the thin-layer Navier-Stokes (TLNS) equations. Due to these forebody 
applications, this code appears to have the capability to address the problem of 
determining the flow field within an inlet of an aircraft at high angles of attack. 
This code also has multigrid capabilities. 

In order to gain some experience in the use and understanding of the code, 
a simple problem is chosen, which is the prediction of the flow about a cone in 
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supersonic flow at an angle of attack. This configuration was chosen because the 
geometry of a cone is simple and a detailed data base is available which includes 
off-body measurements of velocity and flow angles. Although the interest is in 
predicting subsonic, vortical flows, the physics of the vortex development is the 
same for subsonic and supersonic flows. In addition, the use of multigridding to 
accelerate the rate of convergence of the numerical solutions is examined using 
CFL3D (Version 2.1). The solutions obtained are compared with the PARC3D 
code (NASA Lewis Version, Ref. 2) which solves the full Navier-Stokes (FNS) 
equations. Solutions were obtained using a Cray Y-MP computer with compiler 
version 4.0.3. 

This report is divided into several sections. A brief description of the ex- 
periment and data is presented, followed by a discussion of the CFL3D and 
PARC3D codes as well as the computational grids and boundary conditions. 
The results obtained with the CFL3D code are presented for coarse and fine 
grids (one-block and three-block grids respectively), along with a discussion of 
the performance of the multigrid algorithms. The results obtained with the 
PARC3D code using the same grids are compared with the CFL3D results. The 
report ends by presenting some conclusions. 

2.0 Description of the Data 

For comparison with the numerical results, Rainbird obtained useful data 
of surface static pressures [Ref. 3] for a cone with an 18.0 inch base and a 12.5° 
half angle. The off-body data, which includes the Mach number and flow angle 
profiles was obtained by personal communications with Rainbird, are compared 
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with the n um erical results. To the best of the knowledge of the authors the full 
details of this data set have not been published. However, data for the present 
case and several others are contained in Appendix A. Permission to publish 
the data in this report was granted by the Director General of the Institute 
for Aerospace Research, Ottawa, Canada. Rainbird [Ref. 3] indicates that no 
boundary layer trip was used due to the high free-stream turbulence present in 
the wind tunnel. He assumes that transition occurs very close to the cone apex 
(less than 10% of the cone length). 

A 3-hole probe was used to survey the flow field. The probe was kept turned 
into the local mean flow direction and thus enabled measurements of local flow 
angle and pitot pressures. Surveys were conducted at an axial position of 85% 
cone length. The upstream Mach number was 1.8, the angle of attack was 15.75° 
and the flow was turbulent. A diagram illustrating this test is shown in Figure 1. 
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3.0 Numerical Modeling 


3.1 CFL3D Code 

The CFL3D code [Ref. 1] solves the thin-layer Navier-Stokes equations us- 
ing upwind differencing with a total variation diminishing (TVD) scheme and 
employs the Baldwin-Lomax turbulence model [Ref. 4]. The TVD scheme elim- 
inates oscillations due to dispersion errors introduced by the higher order terms 
in the upwind differences by shutting off these higher order terms in regions 
of large flow oscillations. Various options axe available for TVD schemes, flux 
vector- differencing, and upwinding accuracy. The options recommended below 
provided the best results and are used in the solutions presented in this report. 
These options include the min-mod flux limiter for the TVD Scheme, the Roe 
flux difference splitting scheme and third order accurate upwinding. The three- 
factor approximate factorization scheme is used to obtain a block tridiagonal 
system of equations. For the Roe scheme, the equations are diagonalized to ob- 
tain a scalar tridiagonal system of equations that yields a more efficient solver. 
A conservative scheme is employed to transfer information between grid blocks 
[Ref. 5] and multigridding is also available to accelerate the convergence of the 
solution [Ref. 6]. 
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3.2 PARC 3D Code 


The PARC3D code [Ref. 2] solves the full three-dimensional Reynolds- 
averaged Navier-Stokes equations in strong conservation form using the Beam 
and Warming approximate factorization scheme, to obtain a block tridiagonal 
system of equations. Pulliam’s scalar pentadiagonal transformation provides for 
an efficient solver. Like CFL3D, the code uses the Baldwin-Lomax turbulence 
model [Ref. 4]. Its implicit scheme uses central differencing with artificial dissi- 
pation to eliminate oscillations in the solution associated with the use of central 
differences. Trilinear interpolation [Ref. 7] is used to transfer information at the 
grid block interfaces when a multiblock grid is used. 

3.3 Grid 

The effects of grid refinement on the numerical solution were explored using 
two different grids. The grids were algebraically generated with the INGRID3D 
code [Ref. 8], and clustered near the surface using hyperbolic stretching func- 
tions. 

The first grid shown in Figure 2 had dimensions of 29 x 37 x 61 points 
in the streamwise, circumferential, and radial directions, respectively. The first 
axial station is located ahead of the cone due to concerns about locating the 
inflow boundary on the cone. The typical value for y + for the first off-surface 
grid point is approximately 8. The grid is spaced uniformly in the circumferential 
direction at 5° intervals and is packed towards the apex of the cone. 

The second grid, not shown due to resolution problems in reproducing the 
plot, consisted of three blocks with dimensions of 33 x 73 x 73, 33 x 41 x 73, 
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and 33 x 41 x 33 in the streamwise, circumferential, and radial directions for 
Blocks 1, 2 and 3, respectively. The grid block numbering is shown in Figure 
3a. For this grid, the first axial station was located at the cone apex and results 
obtained were not affected by placing the first point at the cone apex. The 
use of three grid blocks was chosen in order to resolve the leeward side vortex 
using Block 1, which was much denser than Blocks 2 or 3. Grid blocks 2 and 
3 have a grid distribution similar to that of the single block grid. A value of 1 
was obtained for the typical y + for the first grid point off the surface in Block 
one; in Block 2 the value was 8. The grid is spaced equally at 1° increments in 
Block 1 and was packed towards the Block 1 interface for Block 2. Note that for 
use in CFL3D, the grids are face-to-face while for PARC3D, the grids overlap 
in order to accommodate the linear interpolation scheme used at the grid block 
interfaces. These interfaces are non-contiguous for both codes. 

3.4 Boundary Conditions 

The upstream and outer radial boundary conditions are held fixed at super- 
sonic free-stream conditions. The flow properties are extrapolated for supersonic 
flow at the downstream exit. At the surface, no-slip, isothermal conditions are 
specified. Slip wall boundary conditions are used along the planes of symmetry. 
These boundary conditions are illustrated in Figure 3a and 3b. 

Isothermal conditions were used because the experiment used a blow-down 
wind tunnel in which the surface temperature variation was less than 5 Fahren- 
heit with a run duration of 20 to 30 seconds. Calculations made using adiabatic 
conditions produced the same Mach number and flow angle profiles as the cal- 
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culations made with isothermal wall conditions. 
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4.0 Results 


In this section several major results will be discussed. The first will deal 
with the effects of grid resolution using CFL3D, and the second the use of multi- 
gridding with CFL3D will then be discussed. Following this dicussion of multi- 
gridding, comparisons of PARC3D solutions with those obtained with CFL3D 
will be presented. 

All of the results reported are derived from PL0T3D format files, which 
use node-centered data. Since the CFL3D code is a finite volume code, the flow 
field is determined at the cell centers of the computational grid and not at the 
grid nodes, as with finite difference codes such as PARC3D. The PL0T3D flow 
and grid files were obtained from the CFL3D code by averaging the adjacent 
cell centers of the grid and using these values at the grid nodes. 

As a preliminary check on the functionality of the CFL3D code and to 
ensure proper problem simulation, a laminar case was studied. The free-steam 
conditions were the same as for the turbulent cases (M^ = 1.8, a = 15.75°). 
The residuals associated with this solution dropped 6 orders of magnitude in 
17,000 iterations and continued to drop. The flow field exhibited a much larger 
cross-flow separation than the turbulent calculations, which is consistent for 
laminar flows. 

Two criteria were used to evaluate the convergence of the turbulent solutions 
which are presented in this report. The first criterion was when the residuals 
reached a constant level, which is typical behavior for turbulent solutions. A 
plot of the density residuals for a constant CFL number of 1 and then a CFL 
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number adjusted from 5 to 1 are shown in Figure 4. The residual values reaching 
a constant level for a CFL number of 1 are shown in Figure 4a. A rapid drop 
and rise in the residuals for the CFL number of 5, seen in Figure 4a, is due to 
the code failing to update local time steps after each iteration until the solution 
is restarted. This behavior is not apparent for a CFL number of 1, (see Figure 
4a), but may be attributable to the solution nearing stability limits at the higher 
CFL number. It did not, however, appear to have an adverse effect on the final 
results. The other criterion was when the change in the boundary layer profiles 
in the vortex region reached a minimum. However, truly steady solutions within 
the vortex region were not obtained. Further discussion of convergence issues 
are presented in the section dealing with multigrid solutions. 

4.1 CFL3D Grid Studies 

In Figure 5, Mach number contours are presented in the plane of symmetry. 
The single and 3-Block grid solutions are very similar. The shock may be slightly 
sharper (closer contours) in the 3-Block grid results due to a few more points 
added in the radial direction. There is a small expansion fan along the leeward 
side of the cone. Mach number contours in the cross-plane are shown in Figures 
6a and 6b. Again, the shock appears slightly sharper in the 3-Block grid solution. 
However, there is a dramatic change in resolution of the leeward side vortex as 
can be seen in the enlarged views of this region shown in Figures 6c and 6d. 
The single block grid does not indicate the vortex presence with the exception 
of a rapidly thickening boundary layer. In contrast, 3-Block grid resolved the 
vortex very well with a distinct region of recirculating flow. 
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The circumferential positions around the cone are defined as 0° on the 
windward side and 180° along the leeward side. Mach number profiles for several 
circumferential stations axe shown Figure 7. These sire taken at 85% of the cone 
length. Along the windward side of the cone, (0° to 90°), the single grid and 3- 
Block grid solutions agree very well with each other and the data. In this region 
the boundary layers are very thin and well-behaved. Along the leeward side of 
the cone the boundary layers begin to thicken and separate in the cross-flow 
direction at approximately 155°. 

As can be seen, the 3-Block grid solutions provide much better agreement 
with the data than the single block grid. In particular, at 170°, the vortex is 
only resolved with the 3-Block grid. It should be noted that the increase in the 
number of grid points solely in this region (block 1) did not improve this result 
very much and this result is not presented. The reduction of the y"*” value of 
the first grid point from 8 to 1 in Block 1 was necessary to provide the results 
shown at 170°. The predicted Mach number profile in the vortex region for 
170° (shown in Figure 8) indicates improved comparison with the data after the 
solution was iterated an additional 10,000 times. The remaining discrepancies 
between the predictions and data may be due to the turbulence model not 
accounting for the vortical flow adequately. This may also be a contributing 
factor to the discrepancies between the predicted and measured Mach number 
profiles at 180°. 

The flow angles are defined in Figure 9, and those predicted with the sin- 
gle and 3-Block grids are compared with data in Figure 10. From the wind- 
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ward plane of symmetry (0°) to 145°, both solutions provided similar results. 
However, in the vortical region from 155°, the 3-Block grid provides improved 
comparisons with the data. 

Surface static pressures at 85% of cone length calculated with both grids 
are compared with data in Figure 11. Very little improvement is shown with 
the 3-Block grid along the windward side of the cone although both grids pro- 
vided good results along the leeward side of the cone. Increasing the number 
of grid points and density of the grid in the radial direction provided for some 
improved shock resolution but offered little improvement in the surface static 
pressure calculations. The good agreement along the leeward side of the cone 
may indicate that inadequate shock resolution on the windward side of the cone 
is the contributing factor in the discrepancies. Some of the discrepancies are due 
to using the difference of two static pressure coefficients in obtaining the coef- 
ficients presented. Decreasing the value for the first off-body grid point in 
Block 2 may improve these results. Although, since the boundary layer profiles 
are in very good agreement with the data, there may be no further improvement. 

Rainbird noted that there was an error in the surface static pressure mea- 
surements due to the windward boundary layer thickness being only twice the 
diameter of the static pressure holes (personal communications). This error di- 
minishes as the boundary layer thickens along the leeward side of the cone. He 
indicated that, the correction to the surface static pressure data was never im- 
plemented because the error was a function of the constantly varying boundary 
layer thickness. This error would account for a small amount of the discrepancy 
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between the calculated and measured surface static pressures along the wind- 
ward side of the cone. Rainbird also indicated that the model alignment error 
was within .1°, therefore, misalignment of the model is probably not an issue. 

One check on the grid dependency of a solution is to compare the velocity 
profiles in unseparated regions with the Law of the Wall. The single and 3- 
Block grid results are shown in Figure 12. Significant improvements were made 
in the comparisons with the fine grid. The discrepancies in the 0° station can be 
attributed to the value of 11 for the y + of the first off-body point, which places 
it out of the viscous sublayer (linear region), making accurate wall shear stress 
calculations impossible. The grid clustering near the wall was not changed in 
Block 2 from the clustering used with the single block grid. 

4.2 CFL3D Multigrid Studies 

Multigridding is a process in which solutions obtained on successively 
coarser grids are used to accelerate the convergence rate for the highest level 
or finest grid. Each successive lower grid has one-half the number of points as 
the next higher level grid. Large scale flow features are developed very rapidly 
with the coarse or lower level grids, while small scale or finer details are resolved 
with the highest level or finest grid because the effectiveness of multigridding 
is problem-dependent; the results reported may not be directly applicable to 
another problem. The results reported in this section are for a single block grid. 

The convergence histories for several multigrid and single grid (non- 
multigrid) schemes are presented in Figure 13, which shows the density residual. 
All of these curves terminate at iterations or cycles where the solution was judged 
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to be converged. Further iterations or cycles, not shown on the plots, did not 
reduce the residual levels further. One criteria used for convergence was when 
the residual histories reached the s am e constant levels. Another criteria used 
to determine convergence was when the Mach number profiles about the cone 
exhibited minimal or no change with additional iterations. This convergence cri- 
teria is illustrated in the selected Mach number profiles shown in Figure 14. As 
can be seen, solutions obtained with all of the schemes used are virtually identi- 
cal with the exception at <j> = 170°. At this location, the possible unsteadiness of 
the vortex may not allow for a truly steady-state solution. Therefore, the point 
where minimum changes in this profile occurred was used as the convergence 
criteria in this region. All solutions were run for 100 iterations in the lammar 
mode prior to running with turbulence. 

The convergence histories of two single grid (non-multigrid schemes) axe 
shown in Figures 13a and 13b. One of these grid schemes used a constant CFL 
number of 1 (Figure 13a). The other one used a CFL number of 5 for 2700 
iterations and then a CFL number of 1 for an additional 1300 iterations (Figure 
13b). As can be seen, the use of a high CFL number for the initial calculations 
reduced the number of iterations required for a converged solution from 8400 
iterations to 4000 iterations. The lower CFL number allows the residuals to 
drop approximately 1.5 orders of magnitude from the level obtained with a CFL 
number of 5. The solution obtained when the residual history became constant 
for a CFL number of 5 is identical to the solution obtained when the CFL number 
was reduced to 1, as shown in Figure 15a. However, the solution obtained after 
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4600 iterations with a constant CFL number of 1 is different for the converged 
solution obtained after 8400 iterations, as shown in Figure 15b. 

The residual histories for two three-level multigrid cycles are shown in Fig- 
ures 13c and 13d. A single three-level V-cycle consists of obtaining solutions on 
two successively coarser grids and then using the corrections obtained from the 
coarse grids to update the solutions on successively finer grids up to the highest 
level. Each three-level W-cycle consists of obtaining solutions on two succes- 
sively coarser grids, using the corrections obtained on the coarse grids to update 
the solution one grid level up, and then return down one grid level. Following 
these coarse grid solutions, the solutions are updated on successively finer grids 
up to the highest level. 

The V-cycle was first run with a CFL number of 5 for 700 iterations and then 
with a CFL number of 1 for an additional 600 iterations. As can be seen from 
Figure 15c, the solution obtained when the residual history became constant for 
a CFL number of 5 is the same as that obtained after reducing the CFL number 
to 1 and iterating until the residuals become constant again. The number of 
multigrid cycles required to obtain a converged answer was 1300, as compared 
with the much larger number of iterations required using the single grids. The 
W-cycle was run with a CFL number of 3 and a converged solution was obtained 
with 700 multigrid cycles. This W-cycle result represents approximately half the 
number of cycles required by the V-cycle to obtain the same level of residual 
drop. 
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4.2.1 Multigrid Results: Single Block Grid 

When examining code performance, several factors which axe shown in Ta- 
ble 1 must be examined. One important factor is the computational speed in 
terms of CPU time per cycle per point. As can be seen in Table 1, the single grid 
(non-multigrid) scheme provides approximately twice the computational speed 
of either multigrid scheme. This difference is due to the additional solutions 
required in each multigrid cycle. However, the actual computational time re- 
quired by these various schemes differed widely. The W-cycle multigrid scheme 
required the least amount of computational time to reach the same level of con- 
vergence as all of the other schemes. In general the multigrid schemes proved 
to be very effective at reducing computational time for this particular problem. 
Part of this effectiveness may be attributable to the fact that the flow had only 
small regions containing three-dimensional effects, specifically the leeward side 
vortex which occupies a very small portion of the flow field. 

4.2.2 Multigrid Results: 3-Block Grid 

The use of multigridding with a three-block grid was also investigated and 
the results are discussed in this section. In order to reduce CPU time the number 
of grid points in the three-block grid was reduced to approximately one-half of 
the original number of points. The computational speed is summarized in Table 
2. As can be seen, the three block grid required significantly more CPU time 
per cycle per point when used in the multigrid mode than the one block grid 
results (shown in Table 1). This increase is attributed to the need to transfer 
information from one block to another in the three-block grid. The convergence 


15 



criteria used for this study was the same as that used for the one-block grid 
multigrid study. Although the multigrid scheme was more costly per cycle, the 
overall time required to obtain a converged solution was reduced significantly 
from the time required for the non-multigrid solution. This result is similar to 
the results obtained using a one-block grid. 
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4.3 CFL3D Comparisons with PARC3D 

One concern in using CFL3D is that the code solves the thin-layer Navier- 
Stokes equations. In this approximation, the derivatives parallel to a surface are 
ignored and therefore regions where there are significant streamwise gradients, 
such as flow reversal, may not be modelled adequately. In order to study the flow 
in this region, the PARC3D code, which solves the full Navier-Stokes equations 
was used to obtain solutions for this cone with the same grids that CFL3D used. 

The computed Mach number profiles obtained with CFL3D and PARC3D 
using the coarse grid (single block grid) are shown with the data in Figure 16. 
The two solutions are identical with the exception of <j> = 180°. It is not clear 
what is causing these discrepancies at this location. In addition, the flow angles 
predicted by the two codes agree well with each other and are shown in Figure 
17. The surface static pressure distributions predicted by CFL3D and PARC3D 
also agree well with each other, as can be seen in Figure 18. For this particular 
case, the thin-layer approximations used in CFL3D do not appear to influence 
the computed results. 

The performance of the two codes for this problem is shown in Table 3. 
The CFL3D code requires approximately 64% more memory than PARC3D. 
The values of the PARC3D storage requirements for thin-layer and full Navier- 
Stokes solutions are instantaneous values displayed during program execution. 
Ideally, the storage requirements are the same. The PARC3D code carries all 
arrays, regardless of solution mode. The CFL3D code was about 50% faster 
than PARC3D in terms of CPU time per iteration per point. However, the 
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actual CPU time required to obtain a converged answer is difficult to state 
since the time required depends on the manner in which the problem is solved. 
For example, running with a large initial CFL number can increase the rate of 
convergence. Another factor affecting the convergence time required, is running 
in the laminar mode for a few hundred iterations, which can reduce the amount of 
computational time significantly. Therefore, several calculations would have to 
be made with each code in order to determine the optimum approach to solving 
this particular problem. However, preliminary comparisons indicate that the 
CFL3D code, when run in the single grid mode, required approximately 42% 
less computational time than PARC3D. 

The performance comparisons for the PARC3D and CFL3D codes using 
the three-block grid are shown in Table 4. The PARC3D memory requirements 
remained about the same, while the CFL3D memory requirement reduced 27% 
as compared to the one-block grid. Overall, the CFL3D code required 15% more 
memory per point than PARC3D using the three-block grid. One advantage of 
PARC3D is that it uses only one grid block in core memory at a time, whereas 
CFL3D keeps all grid blocks in core memory. Therefore, by using additional 
grid blocks, the core memory required by PARC3D can remain constant as the 
number of grid points increases which is not the case with CFL3D. In addition, 
the speed of the codes remained approximately the same as the single block grid 
case and the ratio of total CPU time required for a converged solution using 
CFL3D to PARC3D was similar to the single block grid case. The extra grid 
points in Block 1 of the PARC3D grid are for the required one grid cell overlap 
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which is not needed with the CFL3D code. 


It should be noted that since these comparisons were made, the Cray Y-MP 
compiler was updated to version 5. 0.2.1. For reasons unknown, the speed of the 
PARC3D computations increased approximately 20%. No significant changes in 
the speed of CFL3D were noted. 

Solutions were also obtained with the PARC3D code using the same fine grid 
(3- block) that was used with the CFL3D code. The Mach number and flow angle 
profiles are shown in Figures 19 and 20, respectively. The results obtained with 
the two codes are in excellent agreement with each other. The only significant 
discrepancies occur in the vortex region ( <j> = 155°, 170°). Because the length 
scales used for the Baldwin-Lomax turbulence models in each code are almost 
identical throughout this region, differences in the turbulence models are not 
likely an issue. Differences in the solutions may be attributable to the varying 
amounts of numerical dissipation present in the solutions. The surface static 
pressure distributions obtained with the two codes (shown in Figure 21) axe in 
excellent agreement with each other. The discrepancies between the predictions 
and data have been discussed in a previous section. 

In the process of matching of the turbulence model length scales in the two 
codes, it was found that the search for a length scale was critical to predicting 
the proper location of the vortex. This effect is illustrated in Figure 22 using 
the PARC3D results. The Mach number profile in Figure 22(a) is the result of 
restricting the search for a length scale to the edge of the undisturbed boundary 
layer. This distance happens to correspond to the center of the vortex since the 
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vortex is not much larger than the boundary layer. The profile indicates that 
the predicted vortex position is not the same as the actual position since the 
predicted Mach number profile is different from the experimental profile. The 
Mach number profile shown in Figure 22b is the result of restricting the length 
scale search to the lower edge of the vortex region. This last comparison of 
the predicted profile to the experimental profile improved with this additional 
restriction. The restriction to the lower edge of the vortex eliminated the contri- 
bution of streamwise vorticity to the turbulent viscosity calculations which are 
due to the vortex. Only the transverse component attributable to the attached 
boundary layer was included in the calculation. This result is consistent with 
the original formulation of the Baldwin-Lomax turbulence model. 

The search distance of the turbulent length scale used by the CFL3D code is 
obtained by using the first 64% of the grid points from a surface. The percentage 
of the number of grid points is fixed within the turbulence model subroutines 
and cannot be adjusted by user inputs, as is available in the PARC3D code. 
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5.0 Conclusions 


A major accomplishment of this study was to gain some experience in the 
use of the CFL3D code. The use of block grids and multigridding in analyzing 
the flow about the Rainbird cone has been explored successfully and the appli- 
cation of these techniques to a complicated configuration such as the F-18 inlet, 
forebody, LEX, and wing should be reasonable. 

The grid studies indicate that significant improvements in the prediction 
of details in the flow field can be made by proper selection of grid density and 
proximity to the surface. A major gain in the agreement of the boundary layer 
in the vortical flow region was obtained by placing the first grid point off the 
surface in this region to a distance within a y + value of 1. Outside of this cross- 
flow separation region, improvements were made by increasing the number of 
grid points without decreasing the distance for the first off-body grid point. 

Despite improvements in the boundary layer profiles with increased grid 
resolution, the predictions of the windward surface static pressure distribution 
did not improve. A small portion of the discrepancies may be due to experimen- 
tal errors attributed to the similarity of the windward boundary layer thickness 
to the diameters of the static pressure holes. 

The use of multigridding indicated a significant reduction in the required 
computational time for this problem. However, the effectiveness of multigridding 
is problem- dependent. The use of multigridding with multiple grid blocks also 
showed a significant reduction in CPU time. 

The CFL3D results compared well with PARC3D, indicating that for this 
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problem, the thin-layer approximation is adequate. Further studies with larger 
recirculating flow regions may be necessary to evaluate properly the range in 
which the thin-layer approximation is valid. This study indicates that with 
proper grid resolution, flow field details may be resolved with an algebraic tur- 
bulence model and may not require the use of higher order turbulence models. 
In addition, the use of proper length scales in the algebraic turbulence model is 
critical to obtaining a good prediction of the vortical flow region. 
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Table 1. Performance Summary for Single-Block Grid 
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PARC PERFORMANCE 


CO 

Z 

P 


z 

p 


CO 






X 

CO 





co 

T— 1 

CO 

26 

4 

3.5 



X 

l 











CM 






r-H 





r— | 

CO 




W 

CO 

X 

CO 



o 

X 

p- 

co 

CO 

CO 

23 

4 

i— H 

CO 

2 

j 

b- 

CO 

X 

l 



<1 

s 

X 

Oi 

CM 




rH. 

CM 


-l-= 

a 

• »—* 

o 

a 

C3 

O 

. •— t 

c3 

CD 

a .-t2 


0) 

'-a 

a 

4-3 

'2 

S’ 


CO 

a 

o 

*55 

C 

CD 


CO 

-a 

t-i 

o 


c 

o 

CJ 

CD 

CO 

to 


CO 

(-4 

CD 

"T3 

t 4 co 

O S 

o 

CL hC 

o 

S-< 

p 


Q bO ' — ' 

-a S P 

■n o pu 

OwO 


a3 

-a 


CD 

a 

• ?-H 

Eh 

P 


8 Oh 

O 


O 

P 

tf 

H 

P 

P 

CO 

p 

p 

u 


<q 

co 

lO 


oo 




oo 


c 

■ rH 

O 

cl 

c 

o 
• ^ 

"cl 

, . S-i 

P' a> 

C .t2 


co 

C 

o 

‘co 

i— I 

CD 


Q 


CO 

XJ 

P. 

o 


a 

o 

o 

CD 

CO 


<D 

T3 

2 

4P> 

* 1 ~H 

CJ 

S 5 


CO 

Pi 

CD 

^ „ v 

Pi CO 

a s 

o 

a pc 
o ip 


co 


<D 

bQ 


-a a P 

pop, 

a £ o 


Pi 

Q 

13 

-a 


<v 

a 

• r— « 

H 

P 


8 Cl, 
O 


26 


Table 3. Solution Statistics for Cray Y-MP Computer (Single-Block Grid) 
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Table 4. Solution Statistics for Cray Y-MP Computer (3-Block Grid) 
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Figure 3a. Grid Block Boundaries (Cross-Stream View) 
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Figure 3b. Boundary Conditions (Streamwise View) 
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Figure 6. Mach Number Contours - Cross-Streamwise View 
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Figure 1 1 . Surface Static Pressure Distribution (x/l = 0.85) 
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Figure 15. Mach Number Profiles at (p ~ 170°, (x/l = 0.85) 
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Appendix A 

Test Data for Several Cases 
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Appendix A: Test Data for Several Cases 
= height above surfaces in inches 
= flow angle in degrees relative to cone generator 
= Mach number / edge Mach number 
= static temperature / edge static temperature 
= velocity / edge velocity 

= velocity component parallel to edge velocity / edge velocity 
= velocity component normal to edge velocity / edge velocity 

- OME ) = Mach number * sin (flow angle - edge flow angle) 

- OME ) = Mach number * cos (flow angle - edge flow angle) 

= velocity component parallel to cone generator / edge velocity 

= velocity component normal to cone generator / edge velocity 

= Mach number * sin (flow angle) 

= Mach number * cos (flow angle) 
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boundary layer on windward generator 
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= angle of attack / cone semi angle 
= circumferential angle at which data was taken 
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RUN NO. 4554 MACH NUMBER 1.801 ALPHA= 22.72 AL/THC= 1.818 PHIPP* 0.0 PE /l'0l)=0 . 54 7 1 
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